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ABSTRACT 


The overall problem of stability of unbraced vertical 
cuts and slope with horizontal ground surface is reviewed. 
Three-dimensional aspects of the problem are discussed in 
relation to factors controlling actual mechanism of stability. 

^ = 0 analysis is carried out to ascertain three- 
dimensional stability of vertical cuts and slopes in homo- 
geneous saturated clay for two cases 1 ) strength constant 
with depth 2) Strength increasing with depth. Different 
composite geometrical failure surfaces are assumed and upper 
bound solution is minimised. Dimensionless stability charts 
incorporating three- dimensional aspects are presented and 
their usefulness in design of deep vertical cuts or slopes 
in saturated clays is discussed. 



CHAPTER - I 


IHTRODUCT lOH 

Overall stability problem of temporary unbraced 
vertical cuts and slopes with horizontal ground surface is 
of considerable economic significance as proper evaluation 
of factor of safety could result in appreciable savings in 
earthwork involved. 

Much work has been done in the field of two dimen- 
sional analysis of factor of safety in which the main assump- 
tions are that the failure length of the cut or slope is 
infinite and the soil is homogeneous. Most of these methods 
are based on limit equilibrium techniques. 

However, for the problem like cuts and slopes having 
finite length, the conventional two dimensional analysis will 
give much lower factor of safety than the actual one in which 
the edge effects are also considered. At the edge portion of 
the failure zone, the radius of the shear surface will 
gradually decrease t oward ends , which means that the s hear 
surface is gradually shifted upward at the edges of failure 
and whose overall effect is to increase the factor of 
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safety of the cut or slope . 

Most stability problems, related to man made or 
natural cuts, slopes and embankments are basically three 
dimensional. For cuts and slopes in normally consolidated 
clays and also compacted earth embankments, the short-term 
stability is of great concern. Therefore a = 0 analysis 
has been carried out. As most two-dimensional stability 
analysis in the field is carried out for /S = 0 condition, 
edge effects evaluated here would bring out reserve factor of 
safty available. This is basically a parametric study 
considering various factor influencing three-dimensional 
stability of cuts and slopes. The importance of evalmtion 
of relevant strength characteristics to be used in any sta- 
bility analysis, as brought out by Bjemun (1973) is not 
considered here as the results of two dimensional and three- 
dimensional analyses are compared for same strength parameters. 
In this study the contribution of edge effects to the stability 
is brought out. Two strength variations with depth are 
considered . 

i.e. 1) Strength constant with depth 

2) Strength increasing with depth 



For simplicity circular arc method is used and is 


extended to get three-dimensional factor of safety. The 
basic assumptions governing this method are retained except 
that the shear surface is now no more restricted to an infinitly 
long cylinder but is taken as a finite length of 2L. Three, types 
of shear surfaces are considered. 

1 ) Cylinder with planar edge surface 

2) Cone attached to a cylinder 

3) Ellipsoid attached to a cylinder 

Results obtained from the analyses are represented 
in the form of dimensionless plots. From these plots one 
can easily find out the factor of safety of finite cuts 
or slopes considering the edge effects and the edge portion 
and cylindrical portion at failure in terms of the height 
of the slope. 



CHAPTER - II 


LITERATURE SDRYEY 

Eor any stability analysis, as the factor of safety 
depends very largely on the strength parameters used, they 
should be evaluated correctly. In this thesis this important 
aspect will not be discussed and reference should be made 
to an excellent paper on this topic (Bjerrum 1973) • Many 
procedures, have been developed for plane strain slope stability 
analysis using limiting equilibrium techniques. All these 
procedures differ ^ themselves by the assumption made 

in order to achieve statical determinacy and the particular 
conditions of equilibrium that are satisfied. 

Early invest igat iors of quay wall failures led Swedish 
engineers to the conclusion that an appropriate stability 
analysis could be based on the assumption of circular failure 
surface. The most important advantage of the circular 

failure surf aces is their significant simplification of the 
mechanics of slope stability analysis. 

In 1 918, Eellenius suggested a method, known 
0 = 0 method, for the analysis of short term stahility 
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slopes in homogeneous soil deposits. In this method the 
shear strength is treated entirely due to cohesion. The 
method of analysis is based on the consideration that , the 
sliding surface which requires the greatest amount of cohesion 
for equilibrium of the earth slope is the most dangerous one 
of all. The minimum value of I is found out by trial and 
error. The procedure consists of choosing several circular 
surfaces and picking out the surface which gives the 
minimum value of P. 

j6 = 0 procedure is only applicable to the condition 
in which shear strength is independent of normal stress. 

For soils having both c and 0 fthe assumption of circular 
shear surface is insufficient to achieve statical determinacy. 
To achieve statical determinacy , log - spiral surface of the 
form r = r^ e® 0 choosen, where r is the, radial 
distance of any point on the spiral from the centre point 
of the spiral, r^ is the reference radius and 0 is the 
angle between r and r^. 

The main advantage of the spiral method is that the 
resultant of the normal stress and the frictional component 
of shear stress always passes throxjgh the centre of the 
spiral. So the moment equation will only involve the weight 
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force and the cohesive resistance of the soil. 

Taylor (1937) has suggested a method for stability 
analysis of slopes in soil having both c and j6. The method 
is based on the consideration that, for a circular surface the 
resultant of the normal stress and the frictional component of 
shear resistance (R) will always lie tangent to a circle of 
radius r Sin 0, called the friction circle. But the summation 
of moments about the centre point will involve the normal 
stress distribution and because the unknown number of co- 
efficients required to describe this distribution cannot be 
solved by three equations of equilibrium of statics, the 
problem is highly indeterminate. 

To make the problem statically determinate it is 
assumed that the resultant (R) of the normal stress and the 
frictional part of the shear resistance is concentrated at 
a single paint. This is equivalent to assuming that the 
normal stresses are concentrated at a single point . But for 
a reasonable distribution of normal stress along the shear 
sixrface the resultant normal force must be less than the 
scalar sum of its components. In order to produce the same 
moments as its components, R must lie tangent to a circle of 
greater radius than the friction circle. Taylor (1937) 
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has shown that the underestimate in the radius for the 
circle of tangency by the friction circle procedure is 
equivalent to using the correct radius for some assumed 
stress distribution with some what lower value of 0 . If 
the magnitude and location of the resultant force due to normal 
and frictional shear stress component are assumed to be 
independent of normal stress distribution, then the equivalent 
value of may be calculated from the moment equation. 

The underestimate in by friction circle analysis 
procedure may be represented by the ratio, 


Sin - Sin 

siH-^ 


( 2 . 1 ) 


where, is the mobilized friction angle correspon- 
ding to single concentrated force and is the friction 
angle for some other assumed normal stress distribution. 

The value of K depends on the assumed normal stress distri- 
bution and the angle subtended by the shear surface at the 
centre. 

lamb e and Whitman (1969) have shown that the 
assumption of two concentrated normal forces at each end of 
the shear surface gives an upper bound solution to the 
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factor of safety satisfying the three conditions of equi- 
librium. The solution is determinate since the magnitude of 
two normal forces and the factor of safety are the only 
unknowns to be calculated from the conditions of equilibrium. 

All the methods described earlier are based on 
the equilibrium of the entire mass above the failxxre surface. 
There are many other methods in which the soil mass is devided 
into a number of slices and equilibriimn of each slice is 
considered. The nimnber of equations and unknowns associated 
with complete equilibrium of entire soil mass is summarized 
in table 1 . 

TABLE 1 

EQUATIONS 

n - Moment equilibrium equations for each slice 

n - Vertical force equilibrixmi equations for each 
slice 

n - Horizontal force equilibri-um equations for each 
slice 


3n - 


Total equations 
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UMIOWNS 

1 - Factor of safety 

n — Normal force on the base of each slice 
n - ' location of normal force on the base of each 
slice 

n-1 - Interslice normal forces 

n-1 - Interslice shear forces 

n-1 - Location of inters lice forces 


5n-2 Total unknowns 


The unknown normal force locations on the base of 
each slice may be eliminated from the equation of overall 
moment equilibrium by considering circular shear surface. 

This equation may be written as - 


n 

2 w.r . 
1 


n 

Sin a -E 
1 


s A1 r 


0 


( 2 . 2 ) 


where, w is the weight of individual slice, 
r is the radius of failinre surface, 
a is the base slope of the slice, and 
n is the total number of slices. 
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Substituting s = c’ + (cf' -u) tan 0 ' (2.3) 

n , 

i c’ 1 + (lif-uAl) tan j6 

p = -J_ (2.4) 

i w Sin a 
1 


in which 'Al’ is the length of the base of a particular slice. 

the 

Assumption made in this formulation is that/ normal force ' N' 
and the weight force act through a point at the centre of 
the base of each slice. 

Fellenius (1927) suggested a method with an additional 

assumption that the interslice forces have zero resultant in 

the direction normal to the failure arc of that slice. With 

the 

this assumption /normal force ’ F’ can be represented by 
H = W Cos a . (2.5) 

Substituting this in equation (2.4) 

n 

4 c Al + (w Cos a - u Al) tan 0 

P = ^ (2.6) 

1 w Sin a 
1 


This procedure provides a direct means of calculating 
the factor of safety from the equations of overall moment 
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equilibrium, however neither force nor moment equilibrium 
is satisfied far the individual slice. In this method the 
assumptions regarding side forces involves n-1 equations 
where there are only n-2 unknowns. Hence the system of slice 
is overdeterminate and in general it is not possible to satisfy 
statics. Thus the factor of safety computed by this method 
will have error. 


When the entire slope is in cohesive soil and it is 
appropriate to use undrained strength throughout the slope, 
then the equation (2.6) becomes greatly simplified. 

For a circular shear surface, this equation becomes 


T = 


n. 

^ s hi.1. 
i=1 ^ ^ 

n 

1 w.. Sin 

i=1 ^ ^ 


(2.7) 


If the undrained strength is constant throughout 
the soil, then the numerator is simply s^^. 1 ^ where Ig^ is 
the arc length of the failure surface. Taylor (1948) has 
produced some charts in which s^/i'^H values are plotted 
against depth factor D for different values of slope 
angle, where D is given by 
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D = 


H + h;j 
H 


( 2 . 8 ) 


Here, is the depth of soil below the toe of the slope 
and H is the height of the slope 

Knowing D and slope angle, mobilized shear stress 
required for equilibrium can easily be found out from these 
charts . 

Spencer (196?) has described a procedure for satis- 
fying complete slice equilibrium for a circular shear surface. 
Assinning that the normal forces are located at the centre of 
the base of each slice, spencer achieved statical determinacy 
with the additional assumption that all resultant interslice 
forces are parallel. The (n-1 ) unknown side force incli- 
nations are thereby replaced by a single unknown inclination, 
resulting into 3n unknowns. Although the solution presented 
by Spencer was only directly applicable to a circular shear 
surface, his proced-ure may be extended to slip surface of any 
shape. 

Bishop (1955) has presented a method with as sumptions 
that normal force and weight force act through the mid-point 
of the base of each slice and the interslice forces have zero 
resultant in vertical direction. With these assumptions, 



(2.9) 


factor of safety was given by 
n 

2 cAxj_ + (w^ - u^Axj_) tan ^ 

p _ 

n 

S w. Sin a . 
i=1 ^ 


where , ( a ) 


Cos Cl . i ") + 
^ L 


tan 1^' tan^^i 


j 


(. 2 . 10 ) 


aj_ is the base angle of ith slice, 
is the width of ith slice, 

and Uj^ is the pore pressure at the base of ith slice 

Whitman and Bailey (196?) have shown that numerical 
difficulties may arise in the use of Bishop’s procedure when 
deep failure surfaces are analysed. In this case, the angle a 
may have large negative values such that l/M^a) is zero or 
negative resulting in the normal force at the base of the 
slice being very large or negative, a situation implying 
an unreasonable stress state. Whitman & Bailey have 
indicated that, if the value of 

tan tancc . 

( 1 + F— ^ ) 

for any slice is less than 0.2 the solution should be rejected 
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Merge nstern and Price (1965, 1967) have presented 
a some what different approach to the solution of complete 
slice equilibrium. The assumption made by Morgenstern and 
Price is that the shear and normal forces between the 
slices are related by the expression 

X. = Af (x)E. (2.11 ) 

J J 

where the assumed function f(x) represents the variational 
relationship between the normal and shear component of 
forces (E and X^ respectively) between the slices. The 
parameter is an unknown factor defining the relation 
between X and E in terms of f(x). By assuming f(x) the 
(n-1 ) values of X are replaced by the single unknown. 

Even with this assumption n additional assumptions are 
required for statical determinacy. These additional ass- 
umptions are achieved in Morgenstern and Price procedure 
by assuming that f (x) varies linearly between each interslice 
boundary at which its values are specified. By assuming this 
variation for f(x), the location of normal forces on the 
base of each slice are thus fixed. So the total number of 
unknowns are 3ii. 

Gibson and Morgenstern (1962) suggested a method h 
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for stability amlysis of cuttings in normally consolidated 
clays in which the strength increases with depth. They 
established an expression for factor of safety, based upon 
j6 = 0 analysis, assuming that the soil strength varies 
linearly with depth, as 

where U is stability number, which depends upon the slope 
geometry. 

Many other procedures are available for slope stability 
analysis (J'anbu, Lowe & Karafaith, Janbus Generalized procedure 
of slice, Seed and Sultan etc.). Prom the investigation of 
homogeneous slopes, the conclusions may be drawn with regards 
to the accuracy of the various analysis procedures as 

1) Janbu’s Generalized procedure of slice and 
Spencer's procedure appear to give very similar 
results and are probably the most accurate 
procedures. 

2) The values of factor of safety calculated by 
the Morgenstern and Price procedure would be 
at least as accurate as the values calculated^^ 

■using either Janbu's Generalized procedure of ^ ^ 
slice or Spencer’s proced-ure. 
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3 ) The log-spiral procedure gives results which 

are equally comparable to those by the procedures 
of Janbu and Spencer, however, so far analyses 
using the log-spiral have only been conducted 
for zero pore pressure cases. 

4) The modified Bishop procedure results in slightly 
lower factors of safety than either the generalized 
procedure of slice, Spencer or log-spiral procedures. 
The difference appeals to be less than 4/^ for values 
of ' r^ ’ as large as 0 . 6 . 

5) The values of factor of safety calculated Using 

any of the above procedures with ’r^’ = 0, are 
generally from to 34^ higher than the 

lower bound equilibrium values calculated by 
friction circle method. 

6) The values of the factor of safety calculated 

by the ordinary method of slices are generally 
from li% to 34^ less than the lower bound 
equilibrium values and as much as lower 

than more accurate values calculated by procedure 
like Modified Bishop, These differences correspond 
to analyses without pore pressures. Bor analyses 
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with pore pressure even larger errors will result by using 
the ordinary method of slices. 

In contrast to the voluminous literature on two 
dimensional slope stability, very little work has been 
done regarding three dimensional problems. Sherard et.al. 
( 1963 ) present a method for three-dimensional analysis, 
particularly in case of high dams in narrow valleys, in which 
the length of the dam is devided into a series of segments 
of equal length and the average cross-section of each segment 
is analysed as 2-1) problem. The factor of safety is then 
defined as the ratio of the sum of resisting forces to the 
sum of the driving forces for all segments of dam length. 

The result of this calculation essentially gives a 'weighted' 
average of the stability of the various sections of the 
embankment. 

Xudhbir and Varadarajan (1975) have analysed the Beas 
dam foundation as a three-dimensional problem. Due to the 
nature of shear zones present in the bed rock, the dam along 
its length consisted of blocks of variable resistance, 
lateral resistance between stable and unstable blocks was ' 
evaluated to compute three-dimensional factor of safety. This 
is an effective stress approach and the main difficulty in 
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adopting effective stress approach is the uncertainty about 
the status of lateral stress on edges and the relevant 
strength parameters to be used in the analysis, 

Baligh and Azzouz (1976) extended the concept of 
two dimensional circular arc method to three dimensional 
stability problems, in which some function g(z) is assumed 
to represent the failure surface. The whole length of the 
failure is devided into a number of segments, each of them 
having very small thickness. All the segments are then 
analysed separately as 2-D problem. Then the 3-D factor 
of safety is represented by the ratio of the sum of the 
resisting moments to the sum of the driving moment, for 
all segments in the failure length, 

Bhandari (1970) sioggested a method for stability 
analysis of mud flows considering the edge resistance. In 
this method the normal force on the plane of shearing is 
found out by 

i) the application of the concept of earth pressure 
calculation 

ii) the possible principal shear distribution on 
the side boundary shear. 
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Basal resistance is found out by Morgenstern & Price 
(1965, 1967 ) or by Janbu et.al. (1956) or by the conventional 
method- The total resistance mobilized against sliding is 
expressed in terms of the factor of safety P, which is then 
evaluated by comparing the mobilised resistance with the 
gravity produced shearing forces. It is shown that edge 
effects are of the order of 18?5. 

Statement of Problem : 

In this thesis an approach similar to the one 

be 

adopted by Baligh & Azzouz (1976) will/ adopted to investigate 
edge effects in slope stability analysis. Different types 
of admissible edge surfaces will be investigated to obtain 
an upper bound solution. Cases of constant and variable 
strength with depth will be investigated. 



CHAPiE-E III 


METHOD OP ANALYSIS 


3.1 Two Dimensional Analysis 5 


Before proceeding witli three-dimensional analysis, 
the minimum plane-strain factor of safety, critical centre 
of rotation and the shear surface giving the minimum factor 
of safety ( 2-D ) is to be found out. Por simplicity 
circular arc method has been used to get minimum factor 
of safety and the analysis is carried out for two cases! 

1 ) Constant strength 

2 ) Strength increasing with depth 

Oase-1 ! 


The expression for driving moment per unit 
length is given (Taylor 193?) ass 




VH 


12 


3 


( 1-2 Cot^P + 3 GotX 0otP + 3 Oot aCotX 


- 3 Got a Cot p) 


(3.1.1) 


The expression for resisting moment per unit length is! 

: 11 ^ = 2.a . E^. 0 . 

vhere , a , p , X are as shown in Pig. 3* la 
and ? is the unit weight of soil. 
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H is the height of the slope or vertical cut 
c is the average undrained strength of the soil. 


So the two-dimensional factor of safety can he given hyi 

(3.13) 


^o " 


Case -2 i 

This case has been discussed hy Gibson and 
Morgenstern (1962). The strength at any depth z can be 


represented byj 
kz 


c = 


H 


(3.1.4) 


B’rom Pig. 3.1 , the resisting moment for small length dl 
at depth z on the shear surface is given as; 

d .dl.c 

where , d 1 = R d 9 


dilj = H^dfi 


kz 


H 


(3.1. 4a) 


where 2 = R Sin 9 -R Sin 

Substituting the value of z in eq.uation (3.1.4a). 


dMj 


R^k 

__ 


( Sin9 - Sina^) d9 


^ R^k a 

Hence, III = — • f ^ (SinS - Sinu-) d9 

^ H 
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or, 5.1^ = Cosug - Cosa^ - Sina^ (a^ - )] 

H 

From Fig. 3-1 , it is clear that : 

^2 - = 2a 


and 


R Sin a^ = R Cosa Co^ ^ - R Sin a SinA 

R (0ogx2 “ ) = H CotA 

H = 2 R Sina SinA 
2 

kR 

= — — L H OotA -24a, 'S. (Oosa CosA - Sin a SinA)l 
^ H 


kH^ 


4 Sin^ a Sin^X 


I Cot A + (1-Cota Got A ) aj 

(3. 1.5) 


Ihere as the expression for driving moment remains same 
as given by Equation (3.1. l). Ihe 2-E factor of safety may 
now be given ass 


rCotA + (1-Oota Cot A ) a] 

Sin^a Sin^A (_ 1-2 OotR + 3 C-ot a Cot j3 +30ot a CntA 


or, 


F. 


- 3 


N. 


Cota Cot 
c 

'? z 


3 c 
9 z 


( 3 . 1 . 6 ) 


(3.1.7) 


For normally oonsolidated day, o/ ■:? z is constant for a 
particular soil. So minimization of ET is needed to obtain 
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least value of 


Hepeated trials have been made to obtain least 
value of F^Ci^e. both the cases and these are 

performed by means of a computer programme (Appendix A). 

R represents the radius of the failure surface corresponding 
to ^omin Pi’s senie the position of the centre of 

rotation as shown in the Fig. 3.1a. 

For a cut of 6 meter deep in soil ( •* = 1.77 t/m^) 
having constant strength with depth, the minimum two-dimensional 
factor of safety is obtained as I'omin" '1*56741 

Radius of the critical circle R = 16.25 meter 

n, = 41.0° 

r = 11.52 meter 


For a 9 meter deep cut in soil where strength is 
increasing with depth ( k = 6.85 t/m , = 1.6 t/m^) the 

minimum two-dimensional factor of safety is found as ; 


'omin 


=0.933 

= 586 meters 


=45 

r = 580 meters 


In this case, the failure surface becomes almost a planar 


surface . 
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3.2 I’liree-Dimensional Analysiss 


Vertical cuts 

Getting all these values from two-dimensional analysis, 
one can proceed to perform three-dimensional analysis. 


let the function f = g( 2 ) represents the failure 
surface of revolution at the edge, which is symmetric about 
an axis ( z-axis) through the critical centre and parallel to 
the length of the slope (Fig. 3.1b). Now the three-dimensional 
factor of safety can be written ass 


* 

F 


1^. 2 1q + Resisting moment at the edges 
2 1^ + Driving moment at the edges 


( 3 . 2 . 1 ) 

where, is the resisting moment per unit length corres- 

ponding to the critical circle, is the driving moment 
per unit length corresponding to the critical circle and 
1^ is the half of the length of the cylindrical portion 
of the failure mass. 

To start three-dimensional analysis, it 
is necessary to decide the location of the axis of rotation 
and the type of an edge sin?faoe which will minimise the : 
three— dimensional f actor of safe ty . In this analysis, a 
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rigid body motion of the entire failure mass (cylindrical 
and edge portions) is assumed about a linear axis of 
rotation. 

The critical centre and surface as found bji two- 
dimensional analysis have been assumed applicable for 
three-dimensional analysis also. Two types of edge surfaces, 
namely, conical and ellipsoidal are considered and critical 
values for planar and spherical edge surfaces for a given 
geometry of cut are also investigated. In the first case 
the surface of revolution becomes a cylinder attached to a 
cone and in the second case it becomes a cylinder attached 
to ellipsoid. 

Considering symmetry of the failure zone about 

the 

z = 0 plane ( Fig. 3.1b), let 1^ be the length of/cylindrical 
portion on one side of the failure zone and 1 be the length 
of the edge portion. The length of the cylindrical and edge 
portion ( 1„ and 1 ) are me as u. ed along a line which lies 
on the surface of the slope at shortest distance from the 
axis of rotation of the sliding mass. The position of the 
line along which maximum length of failure occurs will also 
depend on the position of the axis of rotation ('J'| , r). 

For a vertical cut the value of ^ lies between 40° to 45° 
for both the cases ( i.e. constant strength strength 
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increasing with depth). I'rom Fig. 3.1a it is seen that 
point B on the surface of the slope is at shortest distance(r) 
from point 0. 

In the case of a vertical cut, maximum length of 
failure will be measured along a line through B (Fig. 3* la) 
in the direction perpendicular to the plane of the paper 
(i.e. along the line of intersection of top horizontal 
plane and the vertical plane). 

For computation purposes, the edge portion of the 
failure mass is devided into 'n* number of slices, each 
having very small thickness dz, such that n = 1 /dz. 

Now it is needed to represent all the geometrical parameters 
associated with ith slice in terms of known parameters 
1 , R, r and 1-^ - This needs to be done for both types of 
edge surfaces. 

Conical Edge Surface! 

Consider z-axis as the line through the critical 
centre obtained from 2-D analysis and parallel to the length 
of the slope. Z-axis is shown in Fig. 3.2c- Bet C^ be the 
height of the cone. Now in the plane zox* eq.uation of line 
Xq Zq (Fig. 3.2 a) is given by: 


X* 


R 
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x' = R - 


( 3 . 2 . 2 ) 


1 


dx' 


dz 


1 


/ ^ 

=v ■' 72 ' 

^1 


dz 


dz 


or . 


dSi = 


0 


( cj + R-") 


1 


2 ^ 1/2 


(3.2.3) 


0-j_ can be found out from the boundary condition. 

At X* = r, z = 1 (where 1 is the failure length at the 

edge). 

Rrom q.uat ion ( 3 • 2 . 2 ) s 

R1 


r = R - 


0 


1 


IR 


Cl = 


( R-r ) 


(3.2.4) 


In Mg. 3. 2a let z^ be the distance of the middle plane 

of ith slice from x' axis s 

dz 

then, z. = (2 i -- 1) . 

2 

Radius of the circle confining the ith slice can be found out 
frqm equation (3. 2.2) 


i . e « Rj^ ■ S 



(3.2.5) 
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Elliposoidal Edge: 


¥itla orientation of the z and x’ axes same as in 
the previous case, the eq.uation of the ellipse in x’oz plane 
(Eig. 3 *21) is given by: 




( 3 . 2 . 6 ) 


.,2 _ 


R 


Ol^r ( of - 


{ - z- ) 


dz (3.2.7) 


where, z- 


(2i - 1 ) 


Ot can be found out from the condition that , when: 


x' = r, z = 1 


Hence from equation (3.2.6) 


1 
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Or , 


0 


1 


1 E 

(E^-r^)''/^ 


( 3 . 2 . 8 ) 


Let tile radius of the circle confining the ith slice he 
Hence at z = , x' = E^ 

Erom equation (3.2. 6) j 

Ef = h" ( 1 - ) 

R 

or, = Y( gI - zl ) (3.2.9) 

^1 

Eor known values of E and r (obtained from 2-L analysis) 
and for a particular value of 1, 0^ can be found out from 
equations (3.2.4) and (3*2.8) for two types of edge surfaces. 
Knowing E^ ( Equations 3*2.5 and 3.2*9) one can proceed to 
evaluate other parameters , H^, ^ , cc^ (Eig.5.2c)» 

How the equation of the circle confining the ith slice with 
respect to x and y axis iss 

x^ + y^ = E? (3.2.10) 

Equation of line EG ( from Eig. 3.2c) is ? 

X = a ■ 



32 


2 2 *l /2 

- a ) ' (Taking only positive value) 


% = y - 1 = (Rj 


a 


2)1/2 - b 


{3.2. 11) 


where, a = r Oos"^ 
h = r Sinv^ 

From Fig. 3. 2c. 


(3-2. 12a) 
(3.2.12b) 


or 


H. + b 

Sin^B.-ox = ~J=— 

% 

-1 ^ 

/B^ox = Sin ( ) 

% 


Similarly, 


Z%ox = Sin V ( b/Rj_ ) 


a 


1 

2 


-1 + b 1 

Sin ( — ) - Sin (b/Rj_) 


R^ 


(3.2.13) 


Again from triangle 00' Bj_ 


o'Bj_ = oBj_ Sin aj_ 


= \ Sin cc ^ 

Bj_Dj_ = 2 o' B^ =2 R^ Sin® ^ 


^ . = Sin"'* 


H. V -1 

^ ^ = Sin ^ ^ 




2 R^ Sin O^ ^ 


) (3v2.14) 
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Driving and Desisting Moments s 


G ase I { 


Knowing all these values one can represent the 
expression for driving moment for ith slice as : 


'^di 




12 


r 1+3 Cota ^ Cot. -K dz (3-2.15) 


Expression for resiting moment for ith slice, considering 
constant strength iss 




where, c is the strength of the soil. 


( 3 . 2 , 16 ) 


Case II i 

Expression for resisting moment for ith slice 
in the second case (i.e. strength increasing with depth) is s 


^ k 

= -i 

n 


4 Sln^ a^Sin^A^ 


r OotK . + a. (1-Cot A . Cot “.Jhs, 

X 1 1 


(3-2 .17) 

In this case also the expression for driving moment can be 
represented by eq.uation 3*2.15. So finally three- dimensional 


factor of safety can be given by: 

„ M°. 1 + i . M°. 




( 3 . 2118 ) 
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Slopes! 

In case of slopes in parely cohesive soil, the 
failiire length may not always attain maximum value along 
the line of inter section of top horizontal plane and the 
sloping plane, but it may be maximum along any line on the 
surface depending on the position of axis of rotation 
(Fig. 3 *3)* Maximum length of failure will be measured 
along a line which lies on the surface of the slope and at 
shortest distance- from the axis of rotation. 

From Fig. 3*4 a it can be concluded that if o+l3< 90° 
then the point S on the surface of the slope will be at 
shortest distance from ’0’ and consecLuently the failure 
length will attain maximum value along the line through S 
and parallel to the axis of rotation of the failure mass 
(i.e. the projection of the axis of rotation on the inclined 
face of the slope). 

If 90 + P >' 90, then the failure length will 

attain maximum value along the crest of the slope (i.e. line 
throT;igh C and parallel to the axis of rotation), let u be 
the distance from the axis of rotation of the line along 
which maximum length will be attained,- In the first case 
u = OS and in the second case u = 00. Procedure for 






0 •-ritical centre of rotation 

V 

'sX 


n /\ 
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analysing the slopes remains same as for vertical cuts 
except that the expression for different parameters, driving 
and resisting moments are to be modified for both type of 
edge surface ( i.e. conical and ellipsoidal ). 

Oonic al Edge $ 


As before, let is the distance of the middle 

plane of ith slice from the end of the cylindrical portion, so, 

dz 

z. = ( 2 i - 1) — 

^ 2 

' 2 ■ ~ 2 

and, dSj_ = + R 

where, = 1 R / (R-u) (3«2.19) 


Radius of the ith slice is J 


Ri =R 


Rz. 


0- 


Bllipsoidal Edge? 


Here also the expressions for d»j^ , Cj^ and R^ 
will remain same as in the case of vertical cut. 


/ R^ 

dsi - 1 +J-j~ 


4 


he: 


IR 


2 

Cj-Zi 


dz 


^1 = 


2 2 
R^-u^ 


( 5 . 2 . 20 ) 
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and the radius of the ith slices 


^^i = 


2 ^ 

Cl - z. ) 


Knowing and other geometric parameters can be found out. 
Equation of the circle confining the ith slice can be given 
by the equation (3»2.10) i.e.s 

2 , 2 .^2 
X + y = 

Equation of line BO ( Eig. 3 •4b) is given by: 

y = - X tanl^ + ( b + a, tanP* ) 
or y = (a-x)tan'^+b (3*2. 21) 

x=a-(b-y) Got3 


Putting this expression of x in equation (3-2.10) the 
following equation is obtained: 


(b-y)^ Cot^ + a^- 2a (b-y) Got p + y^ = S? 

( 3 . 2 . 22 ) 

This is a quadratic equation in y whose two roots y^i and 
yj^ are as follows: 


A +tB 
Gosec^p 


(3.2.23a) 
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A - V B 

y _ (3.2.23b) 

Coseo p 

where , 

A = Cot[j ( a+b OotP ) (3.2.23c) 

B = R? Gosec^p - (b-Oott^- a)^ (3.2.23d) 


Brom Pig. 3. 4b it is evident that when Rj^ is less than E^, 
i.e. 00 >0S, then both the toe and the heel of failure 
circle will lie on the line BO. Again when R^>Rj^>R^j 
then the toe of the failure circle will lie on line KJ, 
whereas the heel will lie on line 01). Similarly for 
Rj_>R^> R^, the toe of the failure circle will lie on line 
AB and the heel will lie on line OB. Here all the geometric 
parameters are evaluated for the above mentioned cases. 

Oase a s 

Ri < R^ 

Here both the toe(!r^) and heel (P^) of the ith failure 
circle will lie on line BO. 

Prom the Pig. 3.4b it is seen that J 

% = - yi, 

IHiere y^i and y^^ are given by equations (3.2.23a) and (3.2.23b) 
respectively. 

2 fB 

Hence, ~ (3*2. 24a) 

Coec'^p ' 



O! 
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In this case the chord is lying on line BO, 

A i = (3.2.24b) 

Hi 

Chord length 1^ = — = H- Oosec p 

^ ^ sin Ai 



^3 > Rj_ > P4 

In this case as the heel of the ith failure circle is 
lying on line GD, the ordinate of the heel point can be 
given by, y^^ = b. Prom Pig. 3.4bs 

Hj_ = yjc-yji = y^ - 


A + fB 

= — T - 

Oosec^j:S 

( a Ootl:' -b) + E? Cosec^ -(b Ootp- a)^ 
Oosec^ji 

(3.2.25a) 

The equation of line GB can be given by y=b. Substituting 
this in equation ( 3 .2 0 10 ) 5 

X = = (R? - b^)^^^ (3.2.25b) 

Here only the positive value of x is taken since the x- 
co-ordinate of heel of ith failure circle can not be 
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negative. So the co-ordinate of point is 
Substituting in equation (3*2.21 )s 
= - (y(j - "b) Cot p + a 
= a— Hj^ Got p 


So the chord length isJ 


1- P. = 
1 1 


Hence^ 


(^h ^^h 

-1 H. 

Sin . ( — ■— ) 

P.P. 


and 9 


T. P. 

a. = Sin ( — - ) 

2 Ri 


Oase cj 




). 


(3.2.25c) 


(3.2.25d) 


(3.2.25e) 


Rj_ y y R^ 

Here the toe point (I'j^) of ith failure circle will lie on 
line AB, whereas the heel point (Pj_) will lie on line CD. 
bo, as in Oase bJ 

= (R? - b^)l/2 (3.2.26a) 

= b (3.2.26b) 

How the equation of line AB can be given bys 


y = b + H 


(3.2.26c) 


The point of inter section (P^) of line AB and the ith 
failure circle can be given by( 2:^9 y^,). 

Ihere, 

yj, = b + H (3.2.26d) 
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X 

=J- -g? ^ ( 

b+h 

(3.2.26e) 

Chord len 

gths 



liPi 

=/(( x^-x 

J )^ + ( f) 


Hence, 

Sin 

H 

^ I. p. ) — 

11 

(3.2.26f ) 

and 

Sin 

T-P- 

( 

2 R^ 

(3.2.26g) 

hr i ving 

and Resisting Moments 



Knowing all these values we can represent the 
expressions for driving moment of ith slice for both 
the two cases, i.e. Strength constant with depth and strength 
increasing with depth. 

Case Is Constant Strength- 



(3.2.28) 
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Case II! Strengtli Increasing with Depth: 


In this case also the expression for driving moment 
of ith slice can be given by equation (3*2.27). Desisting 
moment, 


I^°i 


k dSj_ 

4 Sih^ aj: Sir^Aj^ 


Got /l^+ a^(l-Cot 
U 



( 3 . 2 . 28 ) 


Finally three-dimensional factor of safety can be given by: 

n 


* 

P = 






1 


c 


+ 



n 

I 

i=1 



(3.2.29) 


where n is the total number of slice at the edge portion. 

Results obtained by using the computer programme (Appendix A) 
* / 

in terms of F / Romin presented in the following 

Chapters. 



OHAPTBR lY 


piss3M‘atior and discussion op HSSULIS 

4*1 Yertical Cuts 

Results obtained from the analysis of vertical cut 
considering three types of failure surfaces namely, cylinder 
with planer edge surface, cylinder attached to a cone and 
cylinder attached to an ellipsoid are represented here 
in the form of dimensionless plots bringing out the effect of 
dimensions of edge portion in comparison with those of cylin^ 
drical part of the failure sirrface. 

In case of conical edge shear surface attached to a 

* / 

cylinder. Pig. 4*1 shows the ratio P /P ^ ^ .- ^ as plotted 

against 1/H lor various values of 1 „/h, for the two cases, i.e. 

strength constant with depth and strength increasing with 
* 

depth. Here P represents the three-dimensional factor 
of safety, Pn ^D in the least two-dimensional factor 

of safety. Pigs. 4.2 and 4.3 represent the results lor 
ellipsoidal edge surface, for the two cases, i.e. strength 
constant with depth and strength increasing with depth respecti- 
vely. Prom these figures following observations may be made t 
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1) Since P A'onia is always greater than unity , 

* 

F always exceeds I’oinin* greater values of l^/E, the 

ratio of P /^omin close to unity indicating that 

the edge portion has very little effect on the overall 
stability of the cut and the problem is simplified to plane- 
strain problem. 

ii) The ratio independent of the 

* 

geometry of the cut and the strength of soil but P and 
^omin iK<3ividually depend on the g:@om@-try and the strength 
of soil. 

iii) In case of cylinder with planen edge surface 
(1/H= 0)' as the additional resiting moment is available 
from the planer edge with no alteration in the driving 
moment, the ratio P /I'ojjjijj becomes maximum for any value 
of 1 q/h (Pigs. 4.1 , 4.2 and 4.3). 

iv) Pailure surface having I^/H greater than 4 can be 
considered to represent plane-strain condition for all 
practical purposes. Por fixed values of I^/E the P reaches 
a minimum at the critical value of 1/H which determines 
the most likely length of failure surface. 

v) It can also be seen that, for the problem of 
vertic al cuts in purely cohes ive homogenous soil, ellipsoidal 
failirre surface gives consistantly lower values of P than 
for conical failure surface a rd is therefore Imost likely 



49 


to represeiat actual failure in the field. Similar conclusion 
has been arrived by Baligh and Azzouz (1976). 

Bor predicting the 3-1 factor of safety of vertical 
cut only ellipsoidal type failure surfaces are considered, 
from Bigs. 4.2 and 4*3 it is seen that for an infinite 
vertical cut it is not possible to predict the actual length 
as the ratio of is continuously decreasing and 

approaches unity with the increasing values of 1 »/h. In ideal 

Ic 

case, the failure length should be infinite ( “jj”" = “ )• 

But in practice failure may have a finite length due to 
non liomogenuity of soil. 

Bor a finite vertical cut also, it is not possible 
* 

to predict B directly from Bigs.' 4.2 and 4-3 as there may 

be many combinations of I^/H and l/H, and it is difficult to 

* 

find which of these will give minimum B . A trial and error 

method is suggested to predict the value of B* for finite vertical 
cut . 

If the total length of the cut is 2L then? 

1^/H + l/H = L/H 

Should alv;ays be satisfied for any combination of 1 q/H and H/H.In a 
plot of 1^/H against 1 /h , equation (4* 1 ) can be represented 
by a 45^ line cutting both 4q/H and 1/h axes at a distaa ce 
I/H from the origin (Big. 4.4 and 4.5 ) 
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As both L and H are known ^ this line can easily be drawn 
in the plot of I^/H against 1/H. Every point on this 
line represents a probable combination of 1-/H and 1/H 

C/ 

satisfying eq,uation 4.1. 

* / 

Contours of equal E /^'omin also drawn in these 
plots. The shape of these curves for the two cases, viz., 
constant strength and strength increasing with depth, are 
represented in Eigs. 4*4 and 4*5 respectively. These curves 
are concave upwards and with increasing values of E 
the curves shift downwords. 

* / 

Contour of minimum E /E^^.,- has also been drawn, 

above which no combination of 1„/H and l/H is possible. 

Erom these figures it is obvious that, the contour of ■ 

* / • 

E which is tangent to the line represented by 

■¥r f 

equation 4"1 will represent the critical value ^ /^omin 

for the finite vertical cut having total length 21. Eor 

a given 45° straight line the appropriate contour of 
* / 

^ zoomin’ tangent to this line, can be interpolated 

with the help of Eigs. 4*2 or 4.3 for the two cases of the 
strengths of the soil. The point of tangency will give 
the critical values of 1/H and 1 q/H. 

■* , 

Eig. 4.6 gives values of E /^omin ^ value of 

L/h = 1.83 using the following shapes of edge surfaces? 
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1 ) Cylinder with planer edge. 

2 ) Cylinder with conical edge surface 

3 ) Cylinder with ellipsoidal edge surface and 

4 ) Cylinder with spherical edge surface. 

It will he seen that the ellipsoidal edge surface gives 

* / 

the minimum value of P /^omin' surfaces considered 

are either linear or convex. Even though one would have 
expected conical surface to give the least value but due 
to different values of 1/H and 1„/H in case of conical and 
ellipsoidal surfaces, the relationship between driving 
and resisting moments gives a smaller value of P /^omin 
for ellipsoidal edge surface. This figure suggests that . 
of the four surfaces considered, ellipsoidal edge surface, 
for the mechanism adopted here, gives the lowest value 

of fV^omin* 

4.2 Slopes ! 

In the case of slopes the 3-1' factor of safety (P*) 
will also depend on the slope angle P . To make a 
comparison, how the slope angle (3 effects the 3-h factor 
of safety (P*) with that obtained for, a vertioal cut, 
a homogeneous clay slope with ^=10°, 30 ° 45° and 60° are 
thoroughly analysed for ellipsoidal edge surfaces only. 
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Por the 10® slope also, plotted against 1/H 

lor different values of 1^/H for two oases i.e. constant 

c 

strength and strength increasing with depth (Fig. 4*7 and 4.8) 
From these figures contours of equal factor of safety- 
are also drawn for the above mentioned cases ( Fig. 4. 9 
and 4.10). The nature of these curves are same as they were 
in case of vertical cuts. 


Using Figs. 4.9 and 4*10 is plotted against 

* / 

F /S’omin'’ 4.11 and 4.12 show these relationships 

for the two cases respectively 

wliere, -f-iS = (Rigs- 4.11 and 4.12) is 

the radius of circle obtained from 2-1) analysis. Results 
for a vertical cut are also indicated in Figs. 4.11 and 4.12. 
The curves L/.iR vs ^ /^omin the intermediate slope 

angle will be within the narrow band provided by l1 ^ 10® and 
P = 90°. From these figures it can be concluded that, slope 
angle has a little influence on the ratio F above 

certain values of L/aR. 


As for example, from Fig. 4.11 and 4.12 it is seen 

F* 


that the difference in the ratio of 


F 


for 1:5 = 10 


omin 


and F =90? for I//4.'5JR = 1 . 5 is approximately 25^ and this 
small difference can be neglected . So for the slopes having 




mm-sm 

‘M'Sm 




>v^ Wm§M 


SLOPE ANGLE=10' 

CONSTANT strength 


2-DirT 









/ r.- : i'j4.Viif, w-l K ' ^J‘ 




'i^.^::;iX:-iih- 

S-X'y. ’’‘l?--*'*^;. < 




&#-fcM''4 






t i^‘ . '•* -X-.i 


.: ; . SLOPE ANGLE -OO^ 

strength;, increasing vVITH DET-'^'^-: 


Q /:-U!mens!Cnai 
’0 ' ;- ’3-'Dj,mensiOfia[ 













_j 

b^B:? 





,/,‘ ' -fcif^f’*’', *h?j /rf^ f'r. i. 1 ''%!■* *' kV ^ ^ '***'V 






m 


!■-■- ;.‘ «■■•(. - •! Wa- 








If!' ^ 




’■‘*•■1 


rfitf#;-®';.' ; ‘ - r-IM-J'-i' . Ll-i.-r 


r^W 


axis of roratio 




jv fjfgtaliofl 
KvfeSfes<J*S4«''*S 













62 


ly- 4&S > 1.5j S' value can be taken corresponding 

to p = 90 ° . 

* - 

In Figs. 4*13 and 4=14 contours of equal F ^ 

have been drawn in plois of 1/H against the slope angle ( p ) 
for the tvro cases, naiaely, ( 1 ) Constant strength, ( 2 ) Strength 
incieasing with depth respectively. From these Figures 
It can be seen that for same values of F "tiie length 

of the edge portion is increasing with the decreasing 
values of slope angle ( p ). However within the range of 
slope angle ( p ) from 45° to 90* this change in the ratio 
of 1/H is not significant but for the values of p less than 
45° the change in the ratio of l/l-I is much more. 

* / 

Obtaining the F from Figs. 4.11 or 4.12 

one can enter in Fig. 4.13 or 4.14 respectively to get the 
value of 1/H for a particular slope angle (p ). Figs. 4. 15 
and 4 . 16 represent the view of the failure surface on a 
plane along which the maximum depth of failure occur, for 
two cases, strength constant with depth and strength 
increasing with depth. These figures are drawn for a 
6 meter deep cut with slope angle varies from i0° to 90 ° 
and the total length of the cut is assumed as 18 meter. 

For example, with p= 60°, the Value of E is found 
as 4 . 5 meter for the case when strength is constant with 
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depth. The value of L/A R is found as 4.05, and from 
Fig. 4.11 the corresponding ^ /^omin found as 1.075* 

With this value of f’*/f'oniin ^ =60° ^ the value of 1/H is 

found as 0.9 from T’ig. 4.13. So the length of the edge 
portion should be 5*4 meter to achieve R /^omin ~ 1.075. 

Similar procedure has been adopted for the second 

case also, i.e. when strength is increasing with depth. 

Pigs. 14.15 and 14*16 bring out a very significant aspect 

of slope failures in cohesive homogeneous clays. It will be 

seen that for flat slopes, (3:5:' 10°, the shape of the failure 

surface is markedly controlled by the edge portions (t — * =1*5 

■^c 

and 1.0 for two cases respectively). Such an observation 
has been reported from the study of mud flows on flat slopes 
by Bhandari, (1970). 

As the slope angle increases beyond 45*^ degrees 
the edge effects are about the same upto p =90°. 



CHAPTER - V 


COHCLUSIOES & RECOmiEHDATIOlB 


Based on the study of effects of shapes of edge 
surfaces in homogenous cohesive clay with streiagth constant 
with depth and increasing with depth, the following main 
conclusions are derived. 

1. For a given cut the 3-1 stability, 

is shown to be controlled by the shape of edge 
surfaces, length of cylindrical and edge portions. 
(Pigs. 4-1, 4.2 and 4-3). 

2. The ellipsoidal edge surface is shown to be the 
most critical one and it is suggested that 

it may represent the field situation (Pig . 4.6). 

3. Contours of equal (Pigs. 4.4 and 4.5) 

are presented and it is shown that all the limiting 
cases lie within the relationship indicated by the 
line AA. It is also seen that in homogeneous 
material the full length of the cut will give 
the minimum value of the ratio ^*/^omin' 
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4. Using the contours of P*/? . a trial and 

error method has been suggested to evaluate 
P*/PojQin ^ finite cut. This is very 
important to practicing engineers in the 
determination of additional stability of cuts 
and slopes due to edge effects. 

5- The effects of slope angle on the value of 

J'*/l’o 2 iin been brought out (Pig. 4.13 and 4.14) 
It is shown that the contribution to stability 
from edge effects is very large for shallow 
slopes and the contribution decreases as the 
slope angle increases. The significance of 
this conclusion for flat slopes is very relavant 
to the study of mud flows on gentle slope. 

The amount of mass involved in failure is also 
shown to depend on the slope angle. 

As this study was aimed at evaluating the effect 
of various parameters on contribution of edge effects to 
stability of cuts and slopes a rigid body mechanism assuming 
the same linear axis of rotation was considered. It is felt 
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that now it is important to consider a more realistic 
mechanism in which the edge portions and the cylindrical 
part are taken to rotate about different centres of rotation. 
Such a study is expected to refine the details presented 
in this thesis, however overall conclusion arrived at are 
not likely to change significantly. 

The other important line of investigation in this 
respect would be an effective stress analysis. 
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